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ABSTRACT 


The purpose of this thesis is to find suitable wavs to design feedback compensators 
for high order systems by using Root-Locus methods. 

As a starting point we will examine a motor amplidyne system and a position con- 
trol system that were previously designed using Bode methods. Then we generalize the 
method and extend it to other systems. 

The final subject of this thesis 1s to design feedback compensators as filters by using 
state feedback coeflicients to define zeros of the filter, then we extend this idea to build 


cascade filters. 
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I. INTRODUCTION 


Modern control svstems may be compensated by placing a suitable filter in either 
the forward (cascade) path or in a feedback path. 

Mainly feedback systems have the advantage that the output follows the command 
more accurately, so the steady state error is less, the bandwidth is greater, and the speed 
of response is faster. Also we can say that the effect of external disturbances and internal 
changes in the parameters or structure of the system 1s significantly reduced. 

‘The conventional and more common way to design a Feedback compensator to 
meet a given set of specifications is the Bode plot method, by determining a suitable gain 
cross-over frequency and I/II curve. 

In this thesis we will try to find some rules to do the same design by using Root- 
Locus methods. To achieve this, first we will examine several designs then trv to put the 
general ideas in a set of rules. 

Simulation studies employed Dynanuc Simulation Language (DSL) and Ewald 


package in the IBM 3033/4381 main frame. 


H. PRELIMINARY EXAMPLES 


A. HIGH GAIN MOTOR AMPLIDYNE SYSTEM 
To study the ideas for feedback compensation first we choose an industrial design 
that exists and is available for comparison. Figure 1 shows the whole compensated 


svstem, where 





Figure |. Original Compensated System 


K, = 14.8 (1) 
T 1209.8 
EE s(s + 57.96)(s + 0.123/9.32) (2) 
0.1375 Us + 62.5(s + 71.65 
H(s) _ 0.1375 (s + 62.5)(s + 714.65) (3) 


(s + 1)(s + 5.26)(s + 625) 


Figure 2 shows the open loop Bode plot for the uncompensated system. The un- 


compensated system is badly unstable. 
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We may use two different kinds of block diagram manipulation to analyze this sys- 


tem. First is the ordinary manipulation shown at Figure 3, where 





Figure 3. . Ordinary Manipulation 


= O > | 
a= + 0.0034)0.066)(s 4- 59.35)(5 + 79.2)(s 6000, 이 © 


Analysis of this transfer function provides the open loop Bode plot, Root-Locus. 
magnified Root-Locus and time response of the system which are given by Figures 4 - 
7. Since there 1s no pole in the right half plane, we can get phase margin and gain margin 


from these plots. 
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The second type of manipulation ts done to get more information using the Root- 
Locus method. The first step for this ts shown at Figure 8. Then with the appropriate 


block diagram reduction the svstem becomes as in Figure 9, where 





Figure 8. First Step For Manipulation 
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Figure 9. Final System For Root-Locus Examination 


17905 


E E 26.98Fj45.19) (6) 


0.00926s^(s -- 62.16)(s + 71.92) 


TELS) = En 1)(s + 5.24)(s + 625.33) 


(7) 

By using the loop transfer function GII(s) to draw the Root-Locus we may examine 
the system for root movements. Since the G(s) has two poles in the right half planc, the 
Bode plot cannot be interpreted from thc point of view of gain margin or phase margin, 
but it shows the stability of tlie system. 

The root locus for the uncompensated system is shown at Figure 10. 

From this point we can proceed to explore the design procedure of the designer. 
First of all, since this design was done using Bode design methods, it will be helpful to 
get the Gl and I/II magnitude curve which is shown at Figure 11 for the original system 
at Figure I 

As we may see from the block diagram there are zeros at the origin. For this system 
we need at least two zeros to keep the svstem type number unchanged. Also the spec- 
ifications show that we want to have a verv high error coefficient for steady state accu- 
racy and a very narrow bandividth. Because of these characteristics, there is a dipole near 
the origin. Ht is possible to see this on the Bode plot of Figure 11. 

The designer could have put the third zero very close to origin but in that case 
building the system might difficult. 

Since other poles are far away compared with these zeros and dipole, these compo- 
nents will decide the system behaviour. 

The Bode plot, Root-Locus and magnified Root-Locus around the origin for the 
loop transfer function GII(s) are given bv Figures 12 - 14. Since we have roots in the 
right half plane, the Bode plot cannot be interpreted for phase and gain margin valucs. 

As a result of this problem we may sav that for high error coeflicient and narrow 


bandwidth svstems we mav need zeros at the origin and a dipole very close to the origin. 
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B. POSITION CONTROL SYSTEM 
As a second example we chose a position control system having an amplidyne and 


a DC armature controlled motor as its power element. Figure 15 shows the original 


compensated system, where 





Figure 15. Original Compensated System 


K, 0.233 (8) 
K, = 0.061 (9) 
12000 

G UY 

S RN 
z 1.885(0.0402s + 1) UA 
G, 0.3335 (12) 
G, — 0.00465? (13) 

E 

6, = 一 一 一 一 一 (14) 


s* + 6.55 + 10 


Figure 16 shows the open loop Bode płot for the uncompensated system. The un- 


compensated system is badly unstable. 
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Figure 16.  Uncompensated System Open Loop Bode Plot 


Again we may use two different kinds of block diagram manipulation to analvze this 


svstem. 


First 1s the ordinary manipulation shown at Figure 17, where 





Figure 17. Ordinary Manipulation 


K, 0.233 


12000(s* -- 6.5s + 10) 
0.005955 + 3.628555f + 247.14455% + 14.455? + 18.85 


~ 


G (16) 
Analysis of this transfer function provides the open loop Bode plot, Root-Locus, 
magnified Root-Locus, and time response of the system which are given by Figures 18 
- 21 
Since there is no root in the right half plane, we can get phase and gain margins from 


this plot. 
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Figure 18. 
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The second tvpe of manipulation is done to get more information using the Root- 


Locus method. The first step for this is shown at Figure 22. 





Figure 22. First Step For Manipulation 


Then with the appropriate block diagram reduction the system becomes as in 


Heure 23, where 





Ar A pon œ , 


Figure 23. — Fifial System For Root-Locus Examination 
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(17) 
.— 0.00125" (s -- 72.39) 


(s + 2.5)(5 + 4) > 


Bv using the loop transfer function GH(s) to draw the Root-Locus we may examine 
the system for root movements. Since G(s) has two roots in the right half plane, the 
Bode plot of this cascade combination cannot be interpreted from the point of view of 
gain margin or phase margin, but it shows the stability of the system. 

The root locus for the uncompensated system 15 shown in Figure 24. From this 
point we can procced to explore the design procedure of the designer. 

First of all, since this design was done using Bode design methods, it will be helpful 
to get the G and 1;H magnitude curve which are shown in Figure 25 for the original 
system in Figure 15. 

As we may see from the block diagram there are three zeros at the origin. As in the 
first example we have to have at least one zero at the origin to keep the system type 
number unchanged. Also since we want to keep the error coefficient unchanged we have 
to have one additional zero at the origm. The reason for the third zero ìs to provide a 
dipole near the origin to make the system have a very high error coefficient for steady 
state accuracy and a very narrow bandwidth. 

For the GH(s), the Root-Locus and a close look around the ongin to the Root- 
Locus are given in Figures 26 and 27. 

Now we mav want to see the effect of removing some of the components from the 
compensator. First, if we remove the third zero the Root-Locus becomes as m 
Figure 28. The svstem is unstable. Now we may want to bend the loci by moving one 
of the poles close to the origin toward left. If we start to move the pole which belongs 
to the dipole at s = -2.5; Figure 29 shows the condition when this pole is at s — -50. 
The system is still unstable. Bv putting this pole further from the zero at s — -72.39, 
we can bend the loci toward the left half plane. We may interpret this as follows : when 
we move the pole to the left we increase the sum of the poles and therefore the sum of 
the roots. Tlns affects the Root-Loci by moving the asymptote centroid to left, so the 
loci from the right half plane cross into the left half plane. Figure 30 shows the Root- 


Loci for this condition. Now we may want to examme the specifications of the svstem 


24 


and compare them with our specifications. Since the dominant roots are far awav from 
the origin the accuracy and damping will decrease, and the bandwidth of the system wiłł 
increase. These do not agree with our specifications. 

After all these trials and sample designs we may summarize the results as a set of 


rules for design of feedback compensation : 


l. To keep the system type number unchanged put at least the same number of zeros 
at the origin as the original system's type number. 


2. To keep both the system type number and the error coeflictent same unchanged, 
the number of zeros at the origin should be N + 1, where N is the type number 
of the uncompensated system. 


3. To have high error coefficient and narrow bandwidth for steady state accuracy in- 
cłude a dipołe close to the origin. Put the zero of the dipołe at the origin for ease 
in physical realization. 

In the next chapter we will apply these rules to several different kinds of problems 


and see whether thev work or not. 
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HI. DESIGNING COMPENSATORS 


A. NON-MINIMUM PHASE SYSTEMS 

Before studying systems that have all poles and zeros in the left half plane, we mav 
want to examine non-minimum phase systems. 1f we have a root or roots in the right 
half plane we are not able to use Bode methods for designing compensators, therefore 
the Root-Locus method will be the onlv useful tool. 


To illustrate we chose the following transfer function, which is the roll mode of a 


vertical take-off aircralt. Figure 31 shows the original system, where 





Figure 31. Original System 


77.187(s + 0.225 + j0.6607) 


= (s — 0.0434)(s + 1.464)(5 + 0.21 + j0.844)(s + 48.35 + jd9.34) Só, 


As we may see from the transfer function, there is a real pole in the right half plane. 
For this reason we can not use Bode methods to compensate the system. The open loop 
Bode plot, Root-Locus and magnified Root-Locus for this system are given by Figures 
32 - 34. 

To find the closed loop roots we make the block diagram manipulations. Then the 


closed loop transfer function becomes : 


77.187(s + 0.225 + j0.6607) 


G1) =". = 0036)(s + 0.21 Z j0.84)(s + 1.474)(s + 48.34 + jd8.72) 


(20) 
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Figure 34. Magnified Root-Locus For Uncompensated System 
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From this point on we may start to think in terms of Root-Locus and loop transfer 
functions. 

Since we will want to have loci from the pole in the right half plane toward the left, 
we have to put either a zero or a pole to the left of this pole. Putting a zero at the origin 
will not help us, because in that case we always will have a root in the right half plane. 
So we should put the compensator poles and zeros somewhere in the left half plane. 
Choosing the specific zero and pole locations is an engineering decision and mostly de- 
pends on the specifications. With this problem we will try to illustrate three different 
conibinations. 


t 


First put a zero close to the origin and a pole far away, both in the left half plane. 


In this case the loop transfer function becomes as in Figure 35, where 





Figure 33. Loop Transfer Function Form 


77.187(s + 0.225 F j0.6607) 


0 (s — 0.036)(s 4- 0.21 F j0.84)(s + 1.474)(s + 48.34 F ¡48.72) = 
K (s +1) 
H, (s) = క (22) 


The Root-Locus and magnified Root-Locus for the loop transfer function are given 
in Figures 36 and 37. As we may see from the plots there is a root-locus segment be- 
tween the pole in the right half plane and the zero at -1. 

To find the value for K,, we should get the tabulated points for the Root-Locus and 
see for which gain all of the roots are in the left half plane. For this particular filter the 
minimum value for K, was 203.7. We selected K, = 500 arbitrarily. For this value, our 


systein becomes as in Figure 38, where 
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77.187(s + 0.225 F j0.6607) 


GE) = 75 = 0:0434)5 + 1.A64)(s + 0:21 F j0-843)(5 + 43:35 F jA9.3) 20 
500(s -+ 1) ( 
H (s) = 450) (24) 





Figure 38. Compensated System 


If we make the block diagram manipulation, we get the system shown in 


Figure 398 vw here 





Figure 39, Resulting Equivalent G 


77.187(s + 0.225Fj0.6617Xs + 50) 


Cea) = GS 30.0269) -- 0235 0.842)(s -- 1.535) s - 49.67) 5 48.41348.63) | 


25) 


The Root-Locus, magnified Root-Locus, open loop Bode plot and time response for 
the compensated svstem are given in Figures 40 - 43. As we may see from the Bode plot 


the svstem is stable but we do not have enough gain. 
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Figure 43. Time Response Without Outside Gain 
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We mav handle this by putting another gain block outside the minor loop. The final 


system then becomes as in Figure 44. 





Figure 44. Final Compensated System 


K mav be selected according to the specifications. We chose K = 100 for illus- 


tration. 


77.187(s + 0.225 F j0.6607) 


యన న యయ య ర నాగా (26) 
DOS SIC cO T3579: S38)05 + 48.35 -F549 34) 
500(s + 1) 
H, (s) = (s + 50) (27) 


The open loop Bode plot and time response for this particular values are given in 
Figures 45 and 46. We may change the loop gain according to the specifications. Effects 


of changing the loop gain will be shown for the third combination. 
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Figure 45. 
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Figure 46. 
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The second tvpe of configuration puts two complex zeros close to the origin and two 
poles far away, all in the left half plane. In this case the loop transfer function becomes 


as in Figure 47, where 





Figure 47. Loop Transfer Function Form 


71.187(s + 0.225 F j0.6607) 


0i) 7 "G7 9036) 4- 0:21 T. 0-83) ^ 1.474)(5 + 48.34 F j48.72) 


(28) 
K.(s +1 F j) 


IO 이 


(29) 


The Root-Locus and magnified Root-Locus for the loop transfer function are given 
in Figures 48 and 49. As we can see from the plots we moved the loci toward the left 


half plane. 
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To find a value for K, , we should get the tabulated points for the Root-Locus and 
see for which gain all of the roots are in the left half plane. For this particular filter the 
minimum value for K, was 2125.8. We selected K2 2500 arbıtrarıly. For this value, our 


system becomes as in Figure 50, where 


77.1870 + 0.225 3. j0.6607) 


G(S) =75 = 0,0334)(s + 1.464)(s + 0,21 £j0843 54835 3;1934 00) 
2500(s + 1 F J) 
PAS = E 20)G 4-30) | 





Figure 50. Compensated System 


If we make the block diagram manipulation, we get the system shown in 


Figure 51, where 





Figure 51. Resulting Equivalent G 


77.1870 + 0.225 X: j0.6617)(s + 20)(s + 30) 


Seal) = S 0.007)(5 + 0.223/0.84)(s + 1.427)(s 4- 24.56.35 - 48.873483). 2 
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The Root-Locus, magnified Root-Locus, open loop Bode plot and time response for 
the compensated svstem are given in Figures 52 - 55. As we may see from the Bode plot 


the svstem is stable but the gain ts too low. 
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Figure 53. Magnified Root-Locus For Final System 
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We may handle this by putting another gain block outside the minor loop. The final 


system then becomes as in Figure 56. 





Figure 56. Final Compensated System 


K may be selected according to the specifications. We chose K = 100 for illus- 


tration. 


77.187(s + 0.225 $ j0.6607) 


06) 7 GG 60334) s 4- 1.46)s + 0,21 F j0.844)(s + 48.35 F jd9.34) 


(33) 
2500(s + 1 F j) 


AS) = (CX 20) 4- 30) 


(34) 
The open loop Bode plot and time response for these particular values are given in 
Figures 57 and 58. We may change the loop gain according to the specifications. Eflects 
of changing the loop gain will be shown with the next combination. 
The third type of combination puts two complex zeros and one real pole close to the 
origin and one pole far away, all in the left half plane. In this case the loop transfer 


function becomes as in Figure 59, where 
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Figure 59. Loop Transfer Function Form 


77.187(s + 0.225 F j0.6607) 


e e E A LL dE 

3 (s — 0.036)(5 + 0.21 F j0.84)(s + 1.473)(s + 48.34 F j48.72) = 
Kkz(s + 2 F jO.5 

1L(s) = Ay + 2 $0.9) (36) 


(s + 1)(s + 50) 


The Root-Locus and magnified Root-Locus for the loop transfer function are given 
in Figures 60 and 61. As we may see from the plots the locus is moved toward the left 


half plane. 
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To find a value for K;, we should get the tabulated points for the Root-Locus and 
see for which gain all of the roots are in the left half plane. For this particular filter the 


minimum value for K; was 81.6. We selected A, = 250 arbitrarily. By this, our system 


becomes as in Figure 62, where 





Figure 62. Compensated System 


77.187(s + 0.225 F 0.6607) 


G) = To = 0.0434 s 1- 1.4645 + 0.21 F 0.843](s + 48.35 F 49.34) 


(37) 


250(s + 2 + 0.57) 


AO 


(38) 
lf we make the block diagram manipulation, we get the system shown in 


Figure 63, where 


77.187(s 4- 0.225 ¥ j0.6617)(s + 1)(s + 50) 


= EZ > LZ. iiie... —- 1 1LnLhL925MqQ'U O0] T, 1 
Ges) — C 0:19) 4- 0.236 F j0.87)(s + 1.52)(s + 49.84)(s + 48.37 4&6) 6 d 





Figure 63. Resulting Equivalent G 


The open loop Bode plot and time response for the compensated svstem are given 
in Figures 64 and 65. As we can see from the Bode plot the system is stable but the gain 


Is too low. 
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Figure 65. 
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We may handle this by putting another gain block outside the minor loop. The final 


system then becomes as In Figure 66. 





Figure 66. Final Compensated Systein 


K may selected according to the specifications. We chose different values for illus- 


tration. 


77.187(s + 0.225 + j0.6607) 


GG) 7 G7 59434 4- 1463) 4- 0.21 X J0.844)s -- 48.35 349.34) 


(40) 


250(s + 2 Y 0.5j) 


AS) = etn cee 7 


(41) 


The open loop Bode plot and time response for K = 100 are given in Figures 67 and 
68. 
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` Figure 68. Final System Loop Gain = 100 


We can change the loop gain according to the specifications. Figures 69 and 70 show 
the time responses for this system when this gain is 500 and 1000. As we may see from 
the time response we can change the amount of overshoot by changing this gain. 

As a result of this problem we can summarize the steps for designing feedback 
compensators for non-minimum phase svstems as following : 

1. Find the closed loop transfer function and get the Root-Locus for this. 
2. Pick the filter poles and zeros according to the Root-Locus rules. 


3. Gct the Root-Locus for the loop transfer function. Tabulate the gain values, and 
select the gain that puts all the roots to the left half plane. 


4. Rearrange the system. Look for the Bode plot and time response. If they are ac- 
ceptable, leave it that way otherwise select another gain outside the minor loop. 
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Figure 69. 
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B. FLEXIBLE ROCKET CONTROL SYSTEM 
As a second example we chose the control system of a flexible rocket which is shown 


in Figure 71. This example was extracted from Ref. |. 


Structure Structure and 
Demod. Filter Servo Rigid Body 





Figure 71. Original Uncompensated System 


G,(s) is the transfer function for the rigid body and is defined as: 


ye) 


a Rİ 일 
333 

SEU 43 

OE 33) e 

G(s) = — > — (44) 

s” + 42.25 + 2750 
(s? + 70s + 4000)(s* + 22s + 12800) i 
Csil) =- ma ~ (45) 
(5 + 305 + 5810)(s” + 30s + 12800) 
And the structure is defined as; 
0.686(s + 53)(s — 53)(s* — 152.25 + 14500)(s? + 153.85 + 14500) 

T(s) = 46) 


(s? -- s + 605)(s? -- 45.5s -- 2660)(s^ -- 2.51s -- 3900)(s^ + 3.995 + 22980) 
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After perfornung all calculations the total system becomes as in Figure 72. Since 


the order of the numerator and denominator are quite large, the transfer function is split 


into three parts to fit on the page. 





Figure 72. Combined Uncompensated System 


(이 = G„ x G Xx GG, where ; 


39134576.25(s -- 35 X j52.7)(s - 11 X: j112.6) 


O7 SETSTJATGAISXAPRDG-33Q 大 0 

ie (s — 21.55 £ j21.79)(s + 72.76 x: j100.24)(s — 83.1 3: 61.82) 48 

027 (s 22.75 x j406.29)(s 4- 1.255 £ j62.44)(s + 1.995 F j151.57)(s — 1.48) (98) 
(s + 11.6 F /24.02) 

Go3 = (49) 


(s + 1.6)(s + 0.5 F j24.6) 


As we may see from the transfer function we have four zeros and a real połe in the 
right half plane. In addition to this most of the poles and zeros which are in the left half 
plane are very close to the imaginary axis. 

The Root-Locus for the open loop uncompensated system is given by Figure 73. 
The system is unstable. Since there is a real pole in the right half plane we can not easily 
interpret the Bode plot for this svstem. 

The next step is to find the loop transfer function as in previous exaniples. By per- 
forming the regular block diagram manipulation the closed loop transfer function be- 


comes G(s) 2 G4(s) x Ga(s) x Ga(s) which is 
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39134576.25(s -- 35 X j52.7)(s - 11 X j112.6)(s + 11.6 F j24) 


_ 5 
AE (5 333)(s + 1.9 F j151.5)(s + 15 F j112.2)(s + 13.8 x j74.5) MY 
人 
a(s) = (s+ 1.5 x j60.4)(s + 19.9 F j52.6)(s + 1 E j24)(s — 0.5 X j4.5) - 
Gas) ^ —— LL (52) 
e (s + 24.8 F j44.17) 


From this point on we may start to think in terms of Root-Locus design. As we saw 
from the loop transfer function we have a pair of complex poles and four zeros in the 
right half plane. In addition to this most of the poles in the left half plane are very close 
to imaginary axis and there are four zeros in the right half plane. So probablv some of 
the roots are going to end up at these zeros. This will give us a gain constraint. 

Since we have poles close to the origin and imaginary axis in the left half plane, 
putting a Zero or zeros at the origin would not help us, because in this case we alwavs 
will have roots im the right half plane. So we decided to use a general type of 
compensator which is two zeros close to the origin and two poles far away. The purpose 
of the zeros was to have loci from the poles in the right half plane terminate on these 
zeros. We made several trials with this type configuration but none of them worked. 
An illustration of this is shown in Figure 74, where G,(s) is the closed loop transfer 


function and 





Figure 74. Illustration of General Configuration 


98.6(s + 0.46 F j4.482) 


Ms) = Te TOD 4- 200) 


(53) 
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The Root-Locus for this case is given in Figure 75. 
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As we may see the system is unstable, because roots in the left half plane were 
moving faster than the roots in the right half plane. 

With this tvpe of problem, since the system has a 12" order numerator and 17" order 
denominator, designing the compensator by just examining the system would be impos- 
sible or require many trials. So we decided to use the function minimization subroutine 
in the DSL package in the mainframe. The parameters used in the subroutine and the 
function minimization program are given in Appendix A and B. 

As the result of the program we found that we have to put the poles of the 
compensator close to the origin and to the left of the zeros which belong to the 
compensator. By putting these zeros and poles in this combination, we make the roots 


in the right half plane move faster. Figure 76 shows the block diagram for the compen- 


sated system, where G(s) is the closed loop transfer function and 





Figure 76. First Resulting Compensator of Function Minimization 


14.72(s 4 0.875 F j2.058) 


20 (s+ 2£j8.343) 


(54) 
The Root-Locus for this compensation is given in Figure 77. By getting the tabu- 
lated root locations, we see that for a Root-Locus gain of 0.068 we get all roots in the 


left half plane. Then our final compensator becomes ; 
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EG 
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7] 


The time response for this case is given in Figure 78. As we may see the system 1s 
stable but the settling time is a little long. According to the specs we may want a faster 


svstem than this. 
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Figure 77. Root-Locus for First FM Result 
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Running the Function Minimization program but changing ¢ and w, to get an out- 


put with the desired settling time, we got the coefficients for this compensator as shown 


in Figure 79. G (s) 1s the closed loop transfer function and 





Figure 79. Second compensator Resulting from Function Minimization. 


8.11(s + 1.1531)(s + 4.878) 


ES > e 00 


(56) 
The Root-Locus for this is given in Figure 80. Again by getting tabulated root lo- 
cations we see that for a Root-Locus gain = 0.124 we get all roots in the left half plane. 


Then our final compensator becomes ; 


S? 4 6.03115 4- 5.625 
Hs) = [22222 (57) 
s + 7.3445 + 45.632 


The time response for this case is shown in Figure 81. As we see the settling time 


and amount of overshoot are within specific limits. 
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Figure 81. 


Time Response for Second FM Result 
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Then our total compensated system becomes as in Figure 82, where G(s) 1s the Open 


loop transfer function and 





Figure 82. Block Diagram of „otal Compensated System 


(s + 1.1531)(s + 4.878) 


తానా... 


(58) 

As a result of this problem we may generalize that, 1f we have components close to 
the imaginary axis in both left and right half plane, we have to use pole-zero couples 
close to each other as a compensator. [n the next chapter building filters by using pole 


placement and state feedback methods will be discussed. 
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IV. DESIGNING FEEDBACK COMPENSATORS USING POLE 
PLACEMENT 


To compensate the systems by state feedback, we may use several computer pro- 
grains to find the feedback coefficients. Sometimes observing or feeding back some of 
the states mav not be possible. For such cases we might have to build an estimator to 
estimate the states. 

The purpose of this chapter is to investigate whether it 1s possible to build a filter 
by using these feedback coefficients as zeros and adding extra poles to the feedback 
path. 

To illustrate this we chose the following plants. 


A. ALL COMPONENTS ARE IN THE LEFT HALF PLANE 


For this case the open loop transfer function of the plant is ; 


1000 


SS s(s + 1)(s + 2)(s + 10) 


(59) 

We want to destgn a feedback compensator (a filter) to stabiltze this system after 
determining the state feedback coefficients. Since the uncompensated system ts fourth 
order, in order to have full state feedback we have to name four roots. 

In this thesis we used the SVS (State Variable Systems) [Ref. 2] program. 

To use the computer program we have to determine n-1 roots (where n is the system 
order) and find the fourth root from the characteristic equation to enter the program. 


Our desired roots are 


$57 —1 72, $5 — —4 (60) 
Then H(s) becomes 
H(s) = s` + 65? + 13s + 20 (61) 
If we get the GH(s) to find the characteristic equation 


1000(s* + 6s* + 13s + 20) 


GH(s) = 
s^ 15 cO 


(62) 
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The characteristic equation then becomes 
కో + 10135? + 603257 + 130203 + 20000 (63) 
And the roots of the characteristic equation are ; 


5, 2 —3.988, 55 4 — —0.997 + 1.996, s, — —1007.02 (64) 


The root locations obtained with this design are within 0.3 % of the desired values 
and the 4" root is at —1007.02 . 

Now we may go to the SVS program to find the state feedback coefficients which 
would place the roots exactly at —4, — 1 xr j2, —1007.02. The state feedback coefficients 
NE ZOt [rom SVS are ; 


H(s) =0.9998s5" + 5.9825? + 12.915 + 19.8 (65) 


These coefficients are negligibly different from the ones we obtained by multiplying 
the desired roots. For accuracy we used the coeflicients that we got from SVS. 
For all practical purposes we may accept the k, which is coefficient of si, as 1. Our 


state feedback compensated svstem then becomes as in ligure 83, where G(s) 1s the 


open loop transfer function and 





Figure 83. State Feedback Compensated System. 


H(s) = s? + 5.9825? 4- 12.915 4- 19.8 (66) 
H(s) in factored form is 
H(s) = (s 4- 0.996 3. j1.993)(s 4- 3.997) (67) 


Root-Loci for the loop transfer function GH(s) are given in Figure 84. 
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To preserve the unity feedback we use block diagram manipulation and the system 


then becomes as in Figure 85. 





Figure 85. The Form of The System With Unity Feedback. 


In Figure 85, G(s) is the open loop transfer function and 


H,(s) = s(s + 2.99 F 1.993) (68) 


To convert the feedback compensator to make 1t more realizable as a filter, we have 
to choose the locations of the extra poles quite far away so they can not aflect the 
transient behaviour of the system. If we put the extra poles into the system form shown 
in Figure 83 we get the required gain easily. The new system then becomes as in 


Figure 86, Where G(s) 1s the open loop transfer function and 


K(s + 3.99)(s + 0.996 F j1.993) 


H,(s) = (s + 300)(s + 400)(s + 500) 


(69) 
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Figure 86. System With Extra Poles. 


To find the value of the K we get the tabulated root locations of the Root-Locus for 


the loop transfer function G/7,(s),which is shown in Figure 87. 
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Root-Locus For Filter Feedback Loop Transfer Function. 


Figure 87. 
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Root-locus gain K to bring the roots where we want them was K = 1.905459 x 10° 
. By using the same manipulation as in Figure 85 if we redraw the block diagram of the 
system) to preserve the unity feedback we get the final compensated system shown in 


Figure 88. 


R + 





Figure 88. Final Compensated System 


G(s) is the open loop transfer function and 


_ 19.8 x 1.905459 x 107 
(i5) 7 $3 300) 4- 4001(5 + 500) ME 


1.905459 x 10’s(s + 2.99 + j1.99) 
HAS == mr eee (71) 
(s + 300)(s + 400)(s + 500) 
The time responses for only state feedback and for the filter compensated case are 
Figures 89 and 90. As we may see from the transient responses since all roots are very 
close to the desired places the time responses are almost exactlv same. The only disad- 


vantage of this configuration may be to burld two identical sets of poles to use in dif- 


ferent places. 


B. NON-MINIMUM PHASE SYSTEMS 

As a second example we chose a non-mininum phase system to check if the idea 
works for both cases. 

For this case the open loop transfer function of the plant is 


m 1000 
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Figure 89. Time Response for Only State Feedback Compensated System 
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Figure 90. Time Response for Filter Compensated System 
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We want to design a feedback compensator (as a filter) to stabilize this svstem after 
determining the state feedback coeflicients. Since the uncompensated system is fourth 
order, in order to have full state feedback we have to name four roots. 

To use the computer program we have to determine n-1 roots (where n is the system 
order) and find the fourth root from the characteristic equation to enter the program. 


Our desired roots are ; 
s2=-1+/l,5=-3 (73) 
Then H(s) becomes 
IIs) = 5° + 5s? + 85 +6 (74) 
If we get the GH(s) to find the characteristic equation 


1000(s + 55? + 8s + 6) 


a s 2915 3 1857 = 405 >) 
The characteristic equation then becomes 
కో + 1021s? + 50185? + 7960s + 6000 (76) 
And the roots of the characteristic equation are 
S, = 2.967, s 3 = —0.988 1.01, s, — —1016 (77) 


The root locations obtained with this design are within 1.1 75 of the desired values 
and the 4^ root is at - 1016. 

Now we may go to the SVS program to find the state feedback coefficients which 
would place the roots exactly at —1 cr Jl, — 3, — 1016. 


The state feedback coeflicients we got from SVS are 
H(s) = 0.99995* + 4.9955* + 7.9525 + 5.956 (78) 


These coefficients again are negligibly different from the H(s) we first accepted but 
we used the ones we found from SVS. 
For all practical purposes we may accept the A, which is coefficient of $ as 1. Our 


state feedback compensated system then becomes as in Figure 91. 
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Figure 91. State Feedback to Compensate The System 


G(s) is the open loop transfer function and 
H(s) = s* 4 4.9955? + 7.9525 + 5.956 (79) 


11(s) in factored form is 


H(s) = (s + 0.995 + j0.998)(s + 2.998) (80) 


The Root-Locus for the loop transfer function GH(S) is given in Figure 92. 
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To preserve unity feedback we may use block diagram manipulation. The next form 
of the system then becomes as in Figure 93 where G(s) is the open loop transfer func- 


tion and 


H (s) = s(s + 2.494 X. 1.316) (81) 





Figure 93. The Form of The System With Unity Feedback. 


To convert the feedback compensator to make it more realizable as a filter, we have 
to choose the locations of the extra poles quite far away so they can not aflect the 
transient behaviour of the system. If we put the extra poles into the system form shown 
in Figure 91 we get the required gain easily. The new system then becomes as in 
Figure 94. 





Figure 94. System With Extra Poles 
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G(s) 1s the open loop transfer function and 


K(s + 2.998)(s + 0.995 ¥ j0.998) 


Hy(s) — (s + 300)(s + 400)(s + 500) 


(82) 


To find the value of the K we get the tabulated root locations of the Root-Locus for 


the loop transfer function GI/,(s),which is shown in Figure 95. 
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Figure 95. 
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The Root-locus gain K to bring the roots where we want them was 
K = 1.1257119 x 10° . By using the same manipulation as in Figure 93 if we redraw the 


block diagram of the system to preserve the unity feedback we get the final compensated 


system shown in Figure 96. 





Figure 96. Final Compensated System 


G(s) is the open loop transfer function and 


5.956 x 1.1257119 x 10” s 


GUE (s + 300)(s + 400)(s + 500) 


1.1257119 x 10’s(s + 2.494 + j1.316) 


(s + 300)(s + 400)(s + 500) (S4) 


H13(s) = 

The time responses for state feedback and for the filter compensated case are given 
Figures 97 and 98. As we may see from the transient responses since all roots are verv 
close to the desired places the time responses are almost exactly same. The only disad- 
vantage of this configuration mav be to build two identical sets of poles to use in dif- 


ferent places. 


C. EXTENDING THE IDEA TO CASCADE COMPENSATION 

The procedure for root placement with the feedback filter used only the loop transfer 
function. We may compensate the system with the same components and same gain bv 
using cascade compensation. If we investigate the systems, the characteristic equations 
are the same with the filter in either the feedback path or in the forward path. Except 


this time instead of the poles, the zeros of the filter are the zeros of the closed loop 
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Figure 97. Time Response For Only State Feedback Compensated System 
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Figure 98. Time Response For Filter Compensated System 
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transfer function but the roots of the closed loop transfer function are in the same places 
as with the feedback compensation scheme. 


To show the result of this idea we chose the same plants. For the first one if we put 


the same filter in the forward path the system becomes as in Figure 99, where 





Figure 99. Cascade Compensated System 


1000 
m aee an 3 
G(s) s(s + 1)(s + 2)(s + 10) (83) 
1.905459 x 10'(s E 1 X j2)(s + 4 
I (s) = 1.205459 x 10 (s + l F j2)(s + 4) (86) 


(s + 300)(s + 400)(s + 500) 


The time response for this system is given in Figure 100. As we may see the time 
responses and oscillatory frequencies are almost identical with the feedback compen- 
sation time response which is shown in Figure 90. In this configuration because of the 
derivatives due to the cascade zeros we have faster rise time. 

We checked the roots of the systems for both the feedback and the cascade config- 


uration. They were identical within 0.37 % difference. 
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Figure 100. 


lime Response for Cascade Compensator 
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The second plant becomes as in Figure 101 1f we put the compensator in the for- 


ward path, where 





Figure 101. Cascade Compensated Second System 


m 1000 
SO AN 두 Sap 
> " R 
EG = .1257119 x 10 (s + Ï F j)(s + 3) (88) 


(s +300)(s + 400)(s + 500) 


The time response for this case is given in Figure 102. Again as we may see the time 
response is essentially the same as that of the feedback configuration which 1s shown in 
Figure 98. Here again we have a faster rise time because of derivatives in the cascade 
filter. 

As we may see from the Figures 99 and 100 we put the filter before the plant in both 
cases and we used the DSL (Dynamic Simulation Language) which behaves the same 
way with the hardware. 

Then we decided to simulate the systems with ALCON (Simulation program for 
personal computers). This program uses arithmethic polynomial calculations to calculate 


the step responses for the svstems. 
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Figure 102. 
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Time Response For Cascade Compensator (Second example) 
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]he resulung time responses from this program for both problems are given in Fig- 
ures 103 and 104. As we see high frequency ripples exist for a short time period which 
we didn't see in the DSL simulation. 

To investigate the reason for this difference we interchanged the blocks in Figures 
99 and 100, that 1s, we fed the error signal into the G block and the output from the G 
block into the H block. For this combination resulting step responses are given in Fig- 
ures 105 and 106 which are identical with the time responses we obtained from ALCON. 

For these DSL simulation cases, in the first case the filter is before plant, the plant 
works as Low-Pass filter and the high frequency ripples are not present, but in the sec- 
ond case the filter follows plant, and the ripples appear in the output. 

We could not decide which one of the simulations is true. According to linear the- 
ory, changing the places of the plants should not aflect to the behaviour of the system. 

To investigate the reason for this difference is left as a recommendation. 

As a result of these examples we may generalize the design steps as follows : 

]. Select N - 1 desired roots. 
Multiply them and find the FI(s). 


Find GHi(s) and get the characteristic equation. 


LE A rJ 


Find the roots of the characteristic equation. If they are acceptably close to the 
desired roots, use H(s) as a state feedback directly, otherwise go to SVS or another 
computer program to find the state feedback coefficients to put the roots exactly 
at the desired places. 


5. Choose the extra poles as far away as possible and put them in the denominator 


of H(s). 


6. Draw the Root-Locus for the loop transfer function Gl1(s) (including the extra 
poles) and get the tabulated data. 


7. Select the appropriate gain to have the system roots in the desired places (generally 
this will be very high gain). 


8. If the purpose is to design a cascade filter, put the resulting filter directly in the 
forward path. Otherwise perform the block diagram manipulation as shown in the 
examples and preserve the unity feedback. 
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ALCON Simulation of First Problem 
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Figure 104. ALCON Simulation of Second Problem. 
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Second DSL Simulation of Second Problem 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In this thesis, the development of a procedure to design the feedback compensators 
by using Root-Locus methods and pole placement methods has been presented. Simu- 
lation results and worked out examples have shown that except for extremely high order 
svstems designing the compensators with the Root-Locus method can be used for anv 
plant. The most important result is to design them using pole placement concepts and 
to be able to use them either in cascade or in the forward path. 

Simulation results have also shown that the compensator can be determined bv a 
Function Minimization subroutine directly from the transfer function of an ideal plant 


based on given specifications. 


B. RECOMMENDATIONS 


]. Further research should investigate methods to satisfv certain specifications while 
designing compensators with the Root-Locus method. 


2. In this thesis, it has been shown that Function Minimization can be used to design 
compensators, but the user should write his own program according to the param- 
eters. A package program mav be written for the main frame to make the procedure 
interactive and faster. 


3. Further research mav investigate the effect of gain variations on the root locations 
while using a pole placement method so one can use directly multiplication of de- 
sired roots as H(s). 
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APPENDIX A. CONSTRAINT PARAMETERS OF FUNCTION 
MINIMIZATION 


Subroutine HOOKE uses the Hooke and Jeeves “direct search” method to find the 
local maximum or minimum of a multi-parameter criterion function, CF. [Ref. 3]. 

The algorithm evaluates CF at a base point, X — (X(1), ..., X(v)), then perturbs each 
parameter in turn, by the amount + S7EP(/), and evaluates CF at each new point. 

If none of these points produces a better value of CF, then the stepsizes are de- 
creased by a factor of BETA, and the process is repeated. 

To use HOOKE, the user must initialize the following arguments and, in the main 
program, CALL  HOOKE(X,STEP,N,ITMAX,CFTOL,ALPFIA BETA CR) 
W,IPRINT,MINMAX) 

All of these arguments must be initialized in MAIN, except for X, CF, Q, QQ, and 
W. Recommended values for ALPIIA and BETA are ALPHA = 2., BETA = 0.5. All 
of the arrays, 1.e., X, STEP, Q, QQ and W, must be declared and dimensioned in MAIN. 

The arguments and their meanings mav be explained as following. 

X = The array of N parameter values. The user must supply the initial guesses, e1- 
ther in the DSL program or in MAIN. 

STEP = An array of dimension N containing the initial step sizes to be used in the 
search. 

N = The number of parameters (a positive integer, at most I5). 

ITMAX = The maximum number of function calls to be performed. 

CFTOL 
terminates ( diflerence between the current value and the previous stage value). 

ALPHA = The factor of (Y - X) which is added to Y to get XNEW; a number 


greater than or equal to 1. 


The error in the criterion function to be reached before the program 


BETA = The stepstze reduction factor; a number between 0 and 1. 

CF = The value of the criterion function. 

Q, QQ, W = Arravs of dimension N, to be used as work space. They must be de- 
clared and dimensioned in the main program. 

IPRINT = An integer flag. O for no intermediate printout, I for intermediate 


printout of X, CF, the number of function evaluation and notification of step-reduction. 
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MINMAX = An Integer flag, 


-1 searches for a minimum, 
mum. 


+l searches for maxi- 
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APPENDIX B. COMPUTER PROGRAM FOR FUNCTION 
MINIMIZATION 


*This program finds the necessary parameters to fit the function to some 
*special transient response. By changing zeta and omgn the standart time 
“response may be generated. 


D COMMON /HANDJ/FLAG,ERRFN,A1,A2,B1,B2 
TITLE Rocket parameters 
* FUNCTION MINIMIZATION 


*Put initial conditions according to the previous assuming parameters. 

INCON A10=2. ,A20=5. ,B10=1. ,B20=100. 

*Define your arrays to simulate your system 

ARRAY P(1),Q(2),AC(1),B(2),D(5),E(5) ,F(1),6(3),H(9),L(TT) KCZYASSE 
UCZ AZ AS AŻ Ge) 

*Put the values for the array coefficients 

TABLE A(1)=333. „,B(1-2)=1. , 333. ;, D(1:5)z1.592: 01534099 50000 MN 
51200000. ,E( 1-5)=1. ,60. 419510. 5383005783650 00 eee 
F(1)=2750. ,6(1-3)=1. ,42.2,2750. ,H(1-9)=1. ,40.652,... 
2949. 22,1048754.7,1.86426E+08,3.9622E+09,1. 17521E+11,... 
2.12662E+12,1: 10017E+14,LEC12110= 55391295 ooo m 
1378867. ,1.83946E+08,5. 2509E+09,3.462441EFI1 2 8516GEFIZNNR 
1.43994E+14 , 1. 06498E+13,-3. 42095E+14 ,K( 1=2)=137375 97972772 
MC1=2)=1. ,30. ,P(1)=25. ,O(1=2)=1. 100, 6 11=1 
S(1-2)=10. ,50. ,T(1-4)=1. ,4. ,8. ,0. 

PARAM R=1. 

“Changing zeta and omgn changes the specifications of your system 

CONST ZETA=0.40, OMGN=6. 0, 1C=0. 0 


INITIAL 
*This part of the program generates the standart time response 
TOTERR=0. 0 


C1=SQRT( 1. O-ZETA* ZETA) 

C2=ONGN* ZETA 

C3=0MGN*C 1 

PHI-ATAN(C1/ZETA) 
* ASSIGN INITIAL GUESS FOR FUNCTION MINIMIZATION 
"This part gives the starting point to the program 

IF (FLAG. GE. O. ) GO TO 10 

A1=A10 

IF (FLAG. GE. 0. ) GO TO 10 

A2=A20 

IF (FLAG. GE. 0. ) GO TO 10 

B1=B10 

IF (FLAG. GE. 0. ) GO TO 10 

B2=B20 

10 CONTINUE 

*Put your variable coefficients here so program can change it 

Y(2)=A1 

Y(3)=A2 

Z(2)=B1 

Z(3)=B2 
*Write normal simulation program for your system 
DERIVATIVE 
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XO-TRNFR(8, 10,0. ,H,L,X2*2. 849) 


X2=TRNFR(0,2,0. ,F,G,X3) 
X3=TRNFR(4,4,0. ,D,E,X4) 
X4=TRNFR(0,1,0. ,A,B,X6) 
X6=15. *ERG 
ERG=ER-X9 
X9=TRNFR(2,2,0. ,Y,Z,X0) 
ER=R-X0 

DYNAMICS 


"Standart time response is generated by this formula 
SID=1.0-(EXP(-C2*TIME)/C1)*SIN(C3*TIME+PHI) 
ERR=X0-STD 
ERRSQ=ERR*ERR 
TOTERR=INTGRL( 0. ,ERRSQ) 

“As a cost function integral of error square is used 

TERMINAL 

“This part finishes the calculation when min value is reached. 

ERRFN = TOTERR 

CONTRL FINTIM=2. ,DELT=0. 01 

END 

STOP 

FORTRAN 

*From now on just declare your step values and amount of variables 

* MAIN PROGRAM FOR FUNCTION MINIMIZATION 

IMPLICIT REAL*8 (A-H,0-2) 
DIMENSION X(4) ,STEP( 4) ,Q( 4) ,QQ(4) ,W( 4) 
C X(1)=kK,X(2)=P INITIAL VALUES IN DSL ROUTINE 
STEP(1)=0. 5DO 
STEP(2)=3. ODO 
STEP(3)=0. 5DO 
STEP(4)=10. 0DO 


N=4 

ITMAX = 100 

CFTOL = 1.0D-6 

ALPHA = 2. ODO 

BETA = 0. 5DO 

IPRINT = O 

MINMAX = -1 

CALL HOOKE(X,STEP,N, ITMAX,CFTOL, ALPHA , BETA, 
* CF,0,00, W, IPRINT,MINMAX) 

STOP 


END 
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